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HIGHLIGHTS 


• Devolatilisation kinetics and analytical Py-GC/MS were conducted on Nigerian teak. 

. Low ash (<1%) and N (<1%) contents and relatively high HHV (19.8 MJ/lcg) for teak was found. 

• Variation in E between 222 and 300 lcj/mol as a function of a. was observed. 

. More than one-third of the chemical constituents were phenolic organic compounds. 
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Devolatilisation kinetics and pyrolytic analyses were carried out on Tectona grandis (teak) using iso- 
conversional methods (Flynn-Wall-Ozawa and Starink) and analytical Py-GC/MS technique respectively. 
Proximate and elemental analyses were also performed and they showed that the ash and C contents and 
the HHV were 0.7%, 49.6% and 19.8 MJ/kg respectively. Biomass sample was subjected to multiple heating 
rates (5-35 K/min) in thermogravimetric experiments and kinetic parameters were evaluated from the 
non-isothermal TGA curves. The activation energy (E) varied between 222 and 300 kj/mol as a function 
of degree of conversion. Similarly, the pre-exponential frequency factor (A) varied between 9.6 x 10 17 
and 9.55 x 10 24 min -1 . Analytical Py-GC/MS showed the presence of C0 2 , acetic acid, furan + 2-butanone, 
levoglucosan, trans-coniferyl alcohol and lignin derivatives. The proportion of phenolic compounds iden¬ 
tified was more than one-third with isoeugenol, acetoguaiacone, and 4-vinylguaiacol showing 
dominance. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Biomass has assumed a pivotal role among several alternative 
energy resources known to humankind. Globally, biomass cultiva¬ 
tion and utilisation for energy and chemical production is being 
pursued vigorously because of its renewable potentials, sustain¬ 
ability prospects and inherent ability to address some environmen¬ 
tal concerns engendered by fossil fuel utilisation. Biomass wastes 
generation often creates disposal and environmental challenges, 
which make them an ideal feedstock for bioenergy production 
(Lasode and Balogun, 2010). Biomass characteristics are 
significantly influenced by climatic conditions, soil types and soil 
fertility treatment in their respective environment. It is therefore 
important to carry out a basic data analysis of regional biomass 
material for specific applications such as Nigerian teak. 

Thermochemical conversion processes such as pyrolysis (pro¬ 
motes fuels and chemical extraction) and gasification (ensures syn¬ 
gas production for power generation) are viable means for the 
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utilisation of biomass residues (Bridgwater, 2003; Sami et al„ 
2001). A major step in these processes is devolatilisation; the liber¬ 
ation of a wide range of volatile species emerging at different rates. 
An investigation into reaction kinetics is very useful for under¬ 
standing the devolatilisation/degradation mechanisms and facili¬ 
tating the design of devolatilisation processes and the sizing of 
reactors. The optimisation of operating conditions is important 
not only for industrial purposes but also for quality control of the 
product. Thermogravimetric analysis (TGA) has become a veritable 
tool for the evaluation of kinetic parameters in thermal degrada¬ 
tion of materials. Prins et al. (2006) and Damartzis et al. (2011) 
carried out studies on weight loss kinetics and kinetic study of 
the pyrolysis of cardoon (a herbaceous plant) via TGA measure¬ 
ments under isothermal and non-isothermal conditions, respec¬ 
tively. The primary advantage of non-isothermal measurement is 
that a single sample is required throughout the entire temperature 
range for the TGA experiment. However, the use of multiple heat¬ 
ing rates has been recommended so as to improve the accuracy of 
non-isothermal techniques (White and Catallo, 2011). 

The devolatilisation characteristics of biomass decomposition 
can be deduced from non-isothermal TGA and the first derivative 
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TGA (DTG) thermograms. These thermograms are also valuable in 
observing the influence of varied heating rates, devolatilisation 
temperature, and volatile matter released during biomass decom¬ 
position. In addition, the application of several model-fitting tech¬ 
niques to TGA data helps in the determination of kinetic 
parameters. A class of such techniques is the model-dependent ap¬ 
proach which relies on a preliminary assumption of the devolatil¬ 
isation model and reaction order (Gai et al., 2013). An accurate 
prediction of the devolatilisation model is crucial and inevitably 
necessitates a prior knowledge of a probable reaction scheme. In 
practice, the activation energy (£) and pre-exponential frequency 
factor (A) are extrapolated from the model of best fit that was sta¬ 
tistically selected among a series of reaction model fitted to non- 
isothermal data (Zhou and Grant, 2004). This method is error 
prone; it produces an unreliable outcome with large discrepancies 
found among the different reaction models (Chen et al., 2006). 
Additionally, the ideal model adopted fails to truly capture the 
events often encountered in heterogeneous reactions (Gai et al., 
2013; Damartzis et al., 2011). 

Suitable alternatives are the iso-conversional methods com¬ 
monly referred to as model-free techniques. They were originally 
conceived for single step reactions (Biagini et al., 2008), which 
would sufficiently capture thermal processes for homogenous 
materials. However, the thermal decomposition of lignocellulosic 
biomass often involves cleavage of chemical bonds in the poly¬ 
meric components following successive and parallel reaction 
routes. For instance, lignin has been reported to decompose over 
a wide temperature range (473-773 K) (Wongsiriamnuay and Tip- 
payawong, 2010). The most reactive of biomass fractions, hemicel- 
lulose starts decomposition at about 523 K, while cellulose begins 
at 570 I< and shows a maximum mass loss rate at about 660 K 
(Gaur and Reed, 1998). Moreover, the interplay of the consequent 
chemical transformation (phase transition from solid to liquid to 
gas) and the physical processes (heat and mass transfer) also intro¬ 
duces complexities into the interpretation of biomass kinetics 
(Damartzis et al., 2011). Apparently, a more robust model other 
than a single step will be required to describe biomass decomposi¬ 
tion. The following improvements were incorporated into iso-con¬ 
versional method for a better simulation of multi-step reactions: 
(1) performance of TGA runs at multiple heating rates and (2) eval¬ 
uation of £ as a function of degree of conversion. The need for an 
assumption on the devolatilisation model is also eliminated (Gai 
et al., 2013). The intricacies of multi-step reaction are thus better 
modelled by iso-conversional techniques providing a valid result 
for the devolatilisation kinetics. 

Lately, the University of florin, Nigeria in an effort to apply the 
principles of sustainable development to the environment insti¬ 
tuted an ambitious programme of planting 100,000 seedlings of 
teak ( Tectona grandis) for a period of five years and it started off 
with 70,383 seedlings on a 57 hectares of land in 2008 (Oloyede 
et al., 2009). Teak is also one of the common tree species used as 
wooden transmission poles and for other structural purposes in 
Nigeria (Babalola, 2008). Teak harvesting and processing evidently 
generate significant quantities of residues that can be utilised for 
thermochemical conversion. 

Pyrolysis is one of the thermochemical conversion processes 
that may be used in harnessing energy potentials of biomass re¬ 
sources. It takes place between 723 and 873 K, essentially for the 
production of condensed liquid (bio-oil) (Azeez et al., 2010). The 
chemical composition of the bio-oil is often of particular interest 
to scientists for its utilisation as biofuels. Analytical pyrolysis gas 
chromatography/mass spectrometry (Py-GC/MS) is a valuable 
technique for prospecting the wide spectrum of volatile species 
that may be found in bio-oil (Greenhalf et al., 2013). Azeez et al. 
(2010) carried out analytical Py-GC/MS on hardwood, that is, iroko 
(Chlophora excelsa ) of African origin. 


Iso-conversional methods have been used extensively for the 
kinetic study of non-woody biomass such as wheat straw (Cai 
and Bi, 2009), bagasse, cashew nut shells, palm stems, coffee husks 
(Wilson et al., 2011), and orange waste (Lopez-Velazquez et al., 
2013). There exists, however, limited literature on kinetics of woo¬ 
dy biomass resource; particularly of tropical origin. The primary 
objective of this study therefore was to conduct reaction kinetics 
using model-free techniques and analytical Py-GC/MS on Teak of 
Nigerian origin in view of the prospective local production of bio¬ 
fuels from its wastes. 

2. Methods 

2.1. Materials 

Teak was obtained from a timber processing plant in December 
2012 at florin, Nigeria (8°29N, 4°32E). Wood chips 

(15 x 10 x 10 mm) from tree lumber, aged between 40 and 
45 years were produced. These were packed in polyethylene bags 
and transported to the Laboratory in the Department of Forest, 
Rangeland and Fire Science, University of Idaho, Moscow, Idaho, 
USA. The biomass samples were then ground in a Thomas Wiley 
Laboratory Mill Model 4, sieved into particles sizes of 0.25- 
0.5 mm and stored in Ziploc bags at room temperature. 

2.2. Methods 

2.2.1. Proximate and elemental analyses, and heating value 

The method according to ASTM Dll02-84 was used for the 
evaluation of the ash content, while the moisture content (MC) 
measurement was performed on a HB 43-S Mettler Toledo mois¬ 
ture analyzer. The fixed carbon (FC) content was calculated by dif¬ 
ference. Elemental analysis was conducted by CE 440 elemental 
analyzer to determine C, H and N contents and O was calculated 
by difference. The Higher heating value (HHV) of biomass samples, 
in triplicate, was determined using a Parr oxygen bomb calorime¬ 
ter model 1341 according to ASTM D5865-04. 

2.2.2. Compositional analysis 

Compositional analysis was conducted using both wet chemis¬ 
try and instrumental analytical methods to obtain the extractive 
and lignin contents according to ASTM D1108-96 and ASTM 
Dll06-96 respectively. 

2.3. Thermogravimetric analysis (TGA) 

Devolatilisation tests were conducted on a thermogravimetric 
analyser (Perkin Elmer TGA-7, Massachusetts, USA) with N 2 purge 
gas (30 mL/min). Samples (2-5 mg) were analysed at linear heat¬ 
ing rates of 5, 15, 25, and 35 K/min. The experiments were carried 
out under non-isothermal conditions; temperature was held for 
1 min at 50 °C and subsequently ramped from 50 to 900 °C. These 
were done in duplicates and good reproducibility was obtained. 

2.4. Kinetic modelling 

The global kinetic equation for the decomposition rate in solid- 
state may be expressed as a product of Arrhenius’ expression and a 
function of the extent of conversion 

^ = Aexp'-£>/(a) (1) 

where A is the pre-exponential frequency factor (min -1 ), £ is 
the activation energy (J/mol), R is the universal gas constant 
(8.314 J/mol K), T is the absolute temperature (K), a. is the extent 
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of conversion and j{oc) is the reaction model which is a function of 
conversion. 

Thermal decomposition of lignocellulosic materials comprise a 
series of complex reactions involving various chemical constitu¬ 
ents occurring simultaneously. This creates a major challenge for 
adopting a simple kinetic modelling. Model-free techniques are 
more robust and are conventionally applied for the determination 
of apparent E and A. These rely on a set of thermogravimetric 
experiments based on multiple heating rates at a particular frac¬ 
tion of conversion. This eliminates the need for a reaction model 
and ‘compensation effect’ in the determination of kinetic parame¬ 
ters. They provide relatively accurate results when compared with 
some other model-fitting techniques (Damartzis et al„ 2011). At 
constant linear heating rate fi = dT/dt, integrating Eq. (1) by separa¬ 
tion of variables gives: 

Fo m =m=A Jl^ )dT (2) 


If T 0 falls well below the temperature where the rate of decom¬ 
position cannot be measured, the lower limit of the integral form 
vanishes. Let x = (- f T ), Eq. (2) becomes 



The term p(x) symbolises the temperature integral and it does 
not have an exact analytical solution. The two possible approxima¬ 
tions adopted for this study were explained in Sections 2.4.1 and 

2.4.2. 

2.4.1. Flynn-Wall-Ozawa's (FWO) method 

The FWO method is one of the notable model-free techniques 
commonly used for the determination of kinetic parameters. The 
FWO method linearises the temperature integral in Eq. (3) using 
Doyle’s empirical approximation (Doyle, 1965) as shown in Eq. (4). 
logp(x) & -2.315 + 0.457x (4) 

Taking logarithm of Eq. (3) and combining with Eq. (4) gives the 
FWO expression in Eq. (5). 

log ' , - 1 ° g ips- 2 - 3 ' 5 -°- 45 4 < 5 > 

By implication a series of TGA experiments at heating rates /?i, 
f} 2 ■ ■ ■ Pi can be conducted and a temperature T) ,■ at a conversion ra¬ 
tio, otj and heating rate, ft, noted. Then a plot of log ft against Tj) , 
for each of j conversion ratio, ai, a 2 • • • otj will give j iso-conversion- 
al lines whose slopes can be calculated from Eq. (6) (Flynn, 1983). 

Slopes 0.4571 (6) 

Thus, the E at each conversion step can be evaluated from the plot of 
log ji against the reciprocal of absolute temperature. An iterative 
procedure was used to obtain the final value of Arrhenius E so as 
to reduce the error introduced by Doyle’s approximation (ASTM 
E1641-07; Flynn, 1983). 

2.4.2. Starink’s method 

Starink (1996) examined two iso-conversional techniques 
(Flynn-Wall-Ozawa and Kissinger-Akhira-Sunose (KAS)) and 
found out that both conform to the expression in Eq. (7). 



where for FWO s = 0, B = 0.457 and for KAS s = 2 and B = 1. Starink 
optimised the values for constants s and B and proposed that 
s=1.8, while B = 1.0037. It was reported that Starink’s method 


was an order of magnitude more accurate than the other two iso- 
conversional methods (Starink, 1996). Hence, Starink’s method 
can be expressed as Eq. (8) 

'"(£)-&-1.0O37T, (8) 

A plot of In (J^) against reciprocal of absolute temperature 
gives a straight line of which the slope corresponds to -1.0037'. 
The E can thus be calculated from the slope of the graph. 

2.4.3. Evaluation of pre-exponential frequency factor (A) 

The evaluation of A was done according to the modified ASTM 
El 641-07 standard. The A was determined as a function of conver¬ 
sion, and based on the assumption that the decomposition obeys 
first-order kinetics. 

2.5. Analytical pyrolysis-gas chromatography/mass spectrometry (Py- 
GC/MS) 

Biomass samples (<0.1 mg) were injected into a quartz capillary 
tube and pyrolysed at 500 °C in an inert atmosphere (He, 
0.125 MPa) using Pyrojector II, (SGE Analytical Science) coupled 
to a GC-MS E i (FOCUS-ISQ, Thermoscientific) instrument. The com¬ 
pounds were separated on a RTx-5 ms (30 m x 0.25 mm 0, Restek) 
column with a temperature programme of 50 to 290 °C (10 min) at 
5 °C/min. The eluted compounds were identified with authentic 
standards, NIST 2008 library matching and by their mass spectra 
(Pillai et al., 2011). 

3. Results and discussion 

3.1. Biomass characterisation 

Biomass characterisation involves proximate, ultimate, compo¬ 
sitional analysis and heating value determination. The proximate 
analysis gives information on the fixed carbon (FC) and the ash 
contents, while the elemental analysis provides carbon (C), hydro¬ 
gen (H), nitrogen (N) and oxygen (O) content determination. The 
compositional analysis concerns the determination of extractives 
(Ext), and lignin (Lig) contents. The basic analysis data from this 
study are compared with results for tropical hardwood from liter¬ 
ature as shown in Table 1. 

Table 1 shows that the ash content for teak (<1%) is much lower 
than that of iroko (3.5%), while the C content for teak (49.6%) is 
higher than the value for iroko (43.9%). It has been demonstrated 
that low ash and high C contents contribute significantly to high 
HHV for lignocellulosic materials (Greenhalf et al., 2013). For iroko 
and teak respectively, the extractive contents are 8.2% and 3.5%, 
while the lignin contents are 26.2% and 36.5%. Telmo and Lousada 
(2011) have established through regression analysis that variation 
in lignin content influences changes in heating value more than 
extractives. The assertions made above are clearly reflected in 
the trend portrayed by the characterisation data in Table 1 with 
teak possessing a comparatively higher HHV than iroko. The values 
of N content for both samples are the same (<1%); an indication 
that low emissions of NO x can be expected if deployed for bioener¬ 
gy applications. 

3.2. Thermogravimetric analysis (TGA) 

Non-isothermal thermogravimetric experiment conducted at 
different heating rates is a viable technique for evaluating kinetic 
parameters of biomass devolatilisation. The application of iso-con¬ 
versional methods, which are model-free, to non-isothermal data 
makes this possible. The accuracy achievable is more reliable when 
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compositional analyses, and heating values for biomass sample. 


Proximate analysis 3 (wt.%) Elemental analysis 3 (wt.%) Compositional analysis' 3 (wt.%) HHV (MJ/kg) 

FC C Ash C H N O c Ext Lig 

Teak 3.8 0.7 49.6 6.3 0.4 43.7 3.5 36.5 19.8 

lroko d 26.1 3.5 43.9 5.3 0.4 46.9 8.2 26.2 16.5 


b Dry-ash-free. 

c Calculated by difference. 

d Azeez et al. (2010), HHV = higher heating value. 


compared with single dynamic measurements. TGA experiments 
for thermal decomposition of teak at varied heating rates CS¬ 
SS K/min) showed a significant mass loss below 398 K. This early 
decomposition stage is as a result of moisture evaporation and pos¬ 
sibly degradation of light organic compounds. The main devolatil¬ 
isation stage in teak, a major step in all biomass thermochemical 
processes, spanned 523-763 K. At this stage, a wide spectrum of 
volatile species is liberated due to the degradation of the essential 
polymeric fractions that make up lignocellulosic biomass. There 
were two distinctive regions. The first region started at less than 
538 I<, rose sharply to maximum peaks for the respective heating 
rates and decayed rapidly till around 678 K. The area below 
538 K has been traditionally assigned to hemicellulose decomposi¬ 
tion because it is the least stable of the polymeric structure of bio¬ 
mass; thereby decomposing earliest (Burhenne et al„ 2013). The 
maximum DTG peaks for teak at differing heating rates are be¬ 
tween 598 and 630 K; these are attributable to cellulose decompo¬ 
sition. Wongsiriamnuay and Tippayawong (2010) observed that 
cellulose decomposes between 598 and 648 K. The second region 
(660-778 K) may be indicative of lignin decomposition because it 
has been noted that lignin degrades over a wider range of temper¬ 
ature (473-773 K) (Wongsiriamnuay and Tippayawong, 2010). Sig¬ 
nificantly, a further thermal cracldng of some resultant organic 
compounds, whose formation may have been precipitated in the 
prior primary reactions, could also be responsible. Biagini et al. 
(2006) in offering an explanation for the weight loss events that 
succeed primary devolatilisation alluded to the degradation of 
chemical species that may have been formed during the preceding 
thermal devolatilisation. 

Table 2 shows the influence of heating rate on devolatilisation 
parameters essentially for the region around the primary peaks. 
As a preliminary step, some of the terminologies adopted from lit¬ 
erature are hereby defined and discussed. T a and T 0 (onset and off¬ 
set temperature respectively) may be defined as intersections of 
the tangent to TGA curve (Gaur and Reed, 1998). These variables 
describe the commencement and ending of the devolatilisation 
process accordingly. T max is the temperature at maximum mass 
loss rate (DTG max ), while a max is the conversion at that tempera¬ 
ture. It was observed that as the heating rate increased, all the 
devolatilisation parameters with the exception of a max , exhibited 
successive shifts toward higher values. For instance at 5 K/min, 


Influence of heating rate on devolatilisation parameters of teak decomposition. 
Heating rate (K/min) T a (K) T 0 (K) T max (K) DTG max (%/min) a max (%) 

5 568 620 603 -7.77 55 

15 584 627 619 -32.5 49 

25 595 639 627 -39.1 49 

35 600 645 630 -56.0 48 


T a (onset temperature) and T„ (offset temperature) are intersections of the tangent 
to TGA curve. T max is the temperature at maximum mass loss rate, a max is the 
conversion at T max . 


T a , To and T max were respectively at 568, 620 and 603 K, while at 
35 K/min they rose to 600, 645 and 630 K. This is in agreement 
with literature where Biagini et al. (2006) noted that the onset 
temperature for rice husk devolatilisation at 5 K/min was 530 K 
and it shifted to 617 K at 100 K/min. This may be linked to heat 
transfer limitations that are usually associated with higher heating 
rates. A temperature gradient is created between the outer and the 
inner layers of the sample shifting toward higher temperature as 
heating rate increases. The increase in value for DTG max is in con¬ 
sonance with the findings by Damartzis et al. (2011) and Kim 
et al. (2010) but contrasts with the finding by Biagini et al. 
(2006) where DTG max tended towards lower values with rise in 
heating rate. Conversion (a max ) was the highest (55%) at 5 K/min 
while it dropped to the lowest value (48%) at 35 K/min. At low 
heating rates, transfer of heat is facilitated and this encourages 
mass diffusion that eventually leads to a greater mass loss. 

3.3. Determination of activation parameters based on iso-conversional 
methods 

The plots of iso-conversional lines in Figs. 1 and 2 were derived 
from the application of Eqs. (5) and (8) to the TGA data. 

The correlation between the logarithmic function containing 
heating rate and 1/Tat a ratios of 0.15 to 0.85 for FWO and Starink 
methods are shown in Figs. 1 and 2 respectively. Table 3 presents 
the values of E calculated from the individual slopes of the a ratio, 
the correlation factor and the relative difference between the two 
methods used. It was observed that the results for the E obtained 
from both methods were in agreement with a difference of 
<0.1%. The relatively low difference between the methods might 
be due to the use of an iterative procedure in the evaluation of E 
for FWO method as earlier mentioned in Section 2.4.1. The high va¬ 
lue of R 2 between 0.15 and 0.65 conversion level revealed that 
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Fig. 2. Determination of E according to Starink method at heating rates of 5,15,25 
and, 35 K/min. 


there was a strong correlation for the most part. Damartzis et al. 
(2011) recorded similarly high values of R 2 for levels between 
0.20 and 0.80 in the application of iso-conversional methods for 
the kinetic study of Cynara cardunculus. The values at low conver¬ 
sion were not included in the calculation of the E because of low 
correlation which can potentially introduce error in the data inter¬ 
pretation. At this stage of conversion it is mainly a physical reac¬ 
tion (drying process) that takes place and the current 
devolatilisation kinetic models are inadequate to define it accu¬ 
rately. Chen et al. (2013) made the same observation about the 
need to adopt a different model; moisture diffusion mechanism 
to obtain a better result for the drying stage. 

For a. between 0.15 and 0.60, the E hovers between 221 and 
260 kj/mol and spikes to about 300 kj/mol at a. level of 0.70. Dam¬ 
artzis et al. (2011) observed a similar trend of variation in values of 
E with increasing conversion. The variation in E with increasing 
conversion may be attributed to biomass heterogeneity in which 
each polymeric fraction possesses a peculiar devolatilisation char¬ 
acteristic. The degree of polymerisation of the reacting pyrolysis 
compounds (secondary reactions), as the thermal decomposition 
progresses, could also influence variation in E. Hence, the devolatil¬ 
isation kinetics can be accurately interpreted by a multi-step 
mechanism (Damartzis et al., 2011). As the reaction progresses, a 
dense carbon-rich structure (char) formed from either a cross-link¬ 
ing of polymer chains, polycondensation or cyclisation reaction of 


Activation energy for conversio 


i range of 0.15-0.85 usir 


; FWO and Starink methods. 


Conversion £ FWO R 2 A (min -1 ) £ Starink R 2 % 

(a) model model Difference 

(kj/mol) (kj/mol) 


15 

20 

25 

35 

40 

45 

50 

55 


70 

75 

80 

85 


221.76 

232.61 
251.75 
259.15 
252.24 
244.13 
242.90 
246.87 
254.67 
263.18 
286.43 
300.58 

237.62 
235.56 


0.993 5.43 x 10 19 221.93 

0.998 2.01 x 10 2 ° 232.51 

0.995 5.43 x 10 21 251.70 

0.998 1.85 x 10 22 259.07 

0.998 2.87 x 10 21 252.24 

0.999 4.39 x 10 2 ° 244.05 

0.999 1.82 x 10 2 ° 243.27 

0.989 5.89 x 10 2 ° 246.77 

0.981 1.86 x 10 21 254.80 

0.955 1.65 x 10 22 263.01 

0.925 1.10 x 10 24 286.27 

0.758 9.55 x 10 24 300.37 

0.683 6.65 x 10 18 237.45 

0.875 9.60 x 10 17 235.35 

0.958 2.29 x 10 19 268.51 


0.993 -0.08 

0.993 0.04 

0.995 0.02 

0.997 0.03 

0.998 0.00 

0.999 0.03 

0.999 -0.15 

0.988 0.04 

0.980 -0.05 

0.952 0.06 

0.920 0.06 

0.746 0.07 

0.664 0.07 

0.865 0.09 

0.955 -0.09 


Table 4 

Analytical Py-GC/MS of teak showing retention time, molecular ion, compounds and% 


Methyl pyruvate 
Furfural 

Furfural alcohol 

2(5)-Furanone 

1,2-Cyclopentanedione 

Phenol 

C4H 6 N 2 0 

3-Methyl-l, 2-cyclopentanedione- 

3- Methyl-phenol 

4- Methyl-phenol 
Guaiacol 

4-Methyl 1,2-benzenediol + unknown 

3- Ethyl-2-hydroxy-2-cyclopenten-l-one 

4- Ethyl-phenol 
Methyl guaiacol 

1,4:3,6 Dianhydro-glucpyranose 
Di-anhydro-hexosan 

4- Vinylphenol,- 
Di-anhydro-hexosan 

5- Hydroxymethyl-2-furancarboxaldehyde 

3- Methoxy-l ,2-benzenediol,- 
Ethyl-guaiacol 

4- Vinylguaiacol 
Syringol 
Eugenol 
Vanillin 
Isoeugenol 
Propylguaiacol 
CioH 12 0 2 
C,oHi 2 0 2 


Levoglucosan 


4-Vinyl-syringol 

Propioguaiacone 

4-Allyl-syringol 

Syringaldehyde 

Cis-coniferyl alcohol 

4-Propenyl-syringol 

Coniferyl aldehyde + acetosyringone 

Trans-coniferyl alcohol 

Hexadecanoic acid 

C, 5 H 14 03 

Methyl-anthroquinone 
Octadecenoic acid 
Octadecenoic acid 
Diplotene 
Sitosterol 
Total unknown 


side chains evolves (Lopez-Velazquez et al., 2013; Zhou et al„ 
2013). Its degradation evidently requires higher energy and this 
may have been responsible for the sudden rise in E. Kim et al. 
(2010) reported a similar finding in their kinetic study of pine 
trees. The A varies linearly with the E. The least and highest values 
of A are 9.60 x 10 17 and 9.55 x lO^min 1 at a. of 0.8 and 0.7 
respectively. 


3.4. Analytical pyrolysis (Py-GC/MS) of teak at 500 °C 

During pyrolytic conversion processes of biomass, volatiles of 
diverse components are liberated from the breakdown of hemicel- 
lulose, cellulose and lignin fractions. Analytical Py-GC/MS 
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technique is used to investigate the components that may be found 
in pyrolytic oil. Table 4 shows the retention time, molecular ion, 
and the % peak area for the various compounds. 

A total of 55 compounds were identified from the Py-GC/MS 
analysis of teak and about 75% were identified. The main com¬ 
pounds were C0 2 (10.4%), acetic acid (6.4%), furan + 2-butanone 
(5.7%), levoglucosan (5.2%), trans-coniferyl alcohol (4.6%). Acetic 
acid and furfural, products of hemicellulose decomposition, were 
among the earliest compounds to emerge (Gao et al., 2013). It is 
well known that levoglucosan is an intermediate product primarily 
from cellulose decomposition. The formation of levoglucosan is 
from the rupture of glycosidic bonds and the rearrangement of glu¬ 
cose (Gaur and Reed, 1998; Gao et al., 2013). The proportion of 
phenolic organic compounds identified was more than one-third; 
however, the major ones are isoeugenol (2.7%), acetoguaiacone 
(2.1%), 4-vinylguaiacol (2.0%), methyl guaiacol (1.6%), guaiacol 
(1.4%), and vanillin (1.2%). These are from the deconstruction of lig¬ 
nin fraction and it is consistent with findings from literature (Pillai 
et al., 2011; Gao et al., 2013). The organic compound that emerged 
at the 50th minute is tentatively identified as diplotene. 

4. Conclusion 

Two different models were employed (Starink and FWO) to 
determine the kinetic data of teak wood and the results were com¬ 
parable. The activation energy and the pre-exponential frequency 
factor varied between approximately 222 and 300kJ/mol, and 
9.6 x 10 17 and 9.55 x 10 24 min 1 as a function of conversion 
respectively. The presence of CO2, acetic acid, furan + 2-butanone, 
levoglucosan, and trans-coniferyl alcohol was detected in signifi¬ 
cant concentration by analytical Py-GC/MS technique. Furthermore, 
isoeugenol, acetoguaiacone, and 4-vinylguaiacol were the main 
phenolic organic compounds identified. Pyrolysed teak thus, serves 
as a good source of extraction of these chemical compounds. 

Acknowledgements 

We would like to acknowledge the technical help of Hui Li with 
the TGA. 


References 

Azeez, AM., Meier, D„ Odermatt, J„ Willner, T„ 2010. Fast pyrolysis of African and 
European lignocellulosic biomasses using Py-GC/MS and fluidized bed reactor. 
Energy Fuels 24, 2078-2085. 

Babalola, F.D., 2008. Overhead wooden transmission pole (WTP) trade in Ibadan, 
Oyo state Nigeria: issues on product quality and standard. Obeche J. 26 (1), 12- 
22 . 

Biagini, E„ Barontini, F., Tognotti, L„ 2006. Devolatilization of biomass fuels and 
biomass components studied by TG/FTIR technique. Ind. Eng. Chem. Res. 45, 
4486-4493. 

Biagini, E„ Fantei, A., Tognotti, L„ 2008. Effect of the heating rate on the 
devolatilization of biomass residues. Thermochim. Acta 472, 55-63. 

Bridgwater, A.V., 2003. Renewable fuels and chemicals by thermal processing of 
biomass. Chem. Eng. J. 91, 87-102. 

Burhenne, L„ Messmer, J., Aicher, T„ Laborie, M„ 2013. The effect of the biomass 
components lignin, cellulose and hemicellulose on TGA and fixed bed pyrolysis. 
J. Anal. Appl. Pyrol. 101,177-184. 

Cai, J.M., Bi, L.S., 2009. Kinetic analysis of wheat straw pyrolysis using iso- 
conversional methods. J. Therm. AnaL. Calorim. 98,325-330. 


Chen, H„ Liu, N„ Fan, W„ 2006. A new method to explain the model dependence of 
apparent activation energy derived from a single nonisothermal dynamic curve. 
Polym. Degrad. Stab. 91,1726-1730. 

Chen, D„ Zheng, Y„ Zhu, X., 2013. In-depth investigation on the pyrolysis kinetics of 
raw biomass. Part I: kinetic analysis for the drying and devolatilization stages. 
Bioresour. Technol. 131, 40-46. 

Damartzis, T.H., Vamvuka, D„ Sfakiotakis, S., Zabaniotou, A, 2011. Thermal 
degradation studies and kinetic modeling of cardoon ( Cynara cardunculus ) 
pyrolysis using thermogravimetric analysis (TGA). Bioresour. Technol. 102, 
6230-6238. 

Doyle, C.D., 1965. Series approximations to the equations of thermogravimetric 
data. Nature 207, 290-291. 

Flynn, J.H., 1983. The iso-conversional method for determination of energy of 
activation at constant heating rates. J. Therm. Anal. 27, 95-102. 

Gai, C„ Dong, Y„ Zhang, T„ 2013. The kinetic analysis of the pyrolysis of agricultural 
residue under non-isothermal conditions. Bioresour. Technol. 127, 298-305. 

Gao, N., Li, A, Quan, C„ Du, L„ Duan, Y„ 2013. TG-FTIR and Py-GC/MS analysis on 
pyrolysis and combustion of pine sawdust. J. AnaL Appl. Pyrol. 100, 26-32. 

Gaur, S„ Reed, T.B., 1998. Thermal Data for Natural Synthetic Fuels. Marcel Dekker, 

Greenhalf, C.E., Nowakowski, D.J., Harms, AB„ Titiloye, J.O., Bridgwater, AV„ 2013. A 
comparative study of straw, perennial grasses and hardwoods in terms of fast 
pyrolysis products. Fuel 108, 216-230. 

Kim, S„ Kim, J., Park, Y„ Park, Y„ 2010. Pyrolysis kinetics and decomposition 
characteristics of pine trees. Bioresour. Technol. 101, 9797-9802. 

Lasode, O.A., Balogun, AO., 2010. Wood wastes generation in Ilorin metropolis: 
problems, management and prospects. In: Proceedings of the 25th International 
Conference on Solid Waste Technology and Management, Philadelphia, PA 
U.S.A 

Lopez-Velazquez, M.A, Santes, V., Balmaseda, J„ Torres-Garcia, E„ 2013. Pyrolysis of 
orange waste: a thermo-kinetic study. J. AnaL Appl. Pyrol. 99, 170-177. 

Oloyede, L, Ayorinde, K.L., Oladele, F.A., 2009. Greening the campus environment: 
The University of Ilorin Experience. 12th General Conference. Association of 
African Universities, Abuja, Nigeria. 

Pillai, K.V., McDonald, A.G., Wagner, F.G., 2011. Developing a model system in vitro 
to understand tracheary element development in Douglas-fir (Pseudotuga 
mensziesii). Maderas Ciencia y tecnologla 13 (1), 3-18. 

Prins, M.J., Ptasinski, K.J., Janssen, F.J.J.G., 2006. Torrefaction of wood part 1: weight 
loss kinetics. J. AnaL Appl. Pyrol. 77 (1), 28-34. 

Sami, M„ Annamalai, K„ Woodridge, M„ 2001. Co-firing of coal and biomass fuel 
blends. Prog. Energy Combust. Sci. 27, 171-214. 

American Standard of Testing and Materials, 2004. ASTM D5865-04,2004. Standard 
Test Method for Gross Calorific Value of Coal and Coke. ASTM International, 
West Conshohocken, PA. 

American Standard of Testing and Materials, 1984. ASTM D1102-84,2007. Standard 
Test Method for Ash in Wood. ASTM International, West Conshohocken, PA 

American Standard of Testing and Materials, 1996. ASTM D1106-96,2007. Standard 
Test Method for Acid-insoluble Lignin in Wood. ASTM International, West 
Conshohocken, PA 

American Standard of Testing and Materials, 1996. ASTM D1108-96,2007. Standard 
Test Method for Dichloromethane (DCM) Solubles in Wood. ASTM International, 
West Conshohocken, PA. 

American Standard of Testing and Materials, 2001. ASTM El 641-07,2007. Standard 
Test Method for Decomposition Kinetics by Thermogravimetry. ASTM 
International, West Conshohocken, PA. 

Starink, M.J., 1996. A new method for the derivation of activation energies from 
experiments performed at constant heating rate. Thermochim. Acta 288, 97- 
104. 

Telmo, C., Lousada, J., 2011. The explained variation by lignin and extractive 
contents on higher heating value of wood. Biomass Bioenerg. 35,1663-1667. 

White, J.E., Catallo, W.J., 2011. Biomass pyrolysis kinetics: a comparative critical 
review with relevant agricultural residue case studies. J. Anal. Appl. Pyrol. 91 
(1), 1-33. 

Wilson, L, Yang, W„ Blasiak, W„ John, G.R., Mhilu, C.F., 2011. Thermal 
characterization of tropical biomass feedstocks. Energy Convers. Manage. 52, 
191-198. 

Wongsiriamnuay, T„ Tippayawong, N„ 2010. Non-isothermal pyrolysis 
characteristics of giant sensitive plants using thermogravimetric analysis. 
Bioresour. Technol. 101, 5638-5644. 

Zhou, D„ Grant, D.J.W., 2004. Model dependence of the activation energy derived 
from nonisothermal kinetic data. J. Phys. Chem. 108,4239-4246. 

Zhou, S., Garcia-Perez, M„ Pecha, B„ McDonald, A.G., Kersten, S„ Westerhof, R„ 2013. 
Secondary vapor phase reactions of lignin derived oligomers obtained by the 
fast pyrolysis of pine wood. Energy Fuels 27,1428-1438. 



